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HE ENGINEERING EXPERIMENT STATION of the Uni- 
versity of Washington, which was established in December, 
1917, performs a three-fold function: 

1) It stimulates, promotes, and coordinates investigations and 
research in all fields of engineering. 

2) It publishes results of significant research achievements. 

3) It provides opportunities through Graduate Research Assist- 
antships for engineering students to extend their profes- 
sional education while pursuing a course of study leading to 
the M.S. or Ph.D. degree. 

For administrative purposes the work of the Station is organ- 
ized under six divisions. All research is carried on in the various 
academic departments of the College of Engineering under the 
supervision of the department faculties. 

The control of the Station is vested in the Station Board con- 
sisting of the Dean of the College of Engineering as chairman, 
the Director of the Station, and the executive officers of the aca- 
demic departments. 

The Station offers a substantial number of Research Assistant- 
ships to highly qualified graduate students who are assigned to 
the academic departments. Results of major investigations are 
published by the Station in the form of bulletins and reports. 
Current research findings, as well as listings of sponsored research 
projects, appear in the quarterly journal, The Trend in Engineer- 
ing at the University of Washington. Reprints of articles by 
members of the engineering faculty and graduate students pub- 
lished in recognized technical journals are also issued by the 
Station. 

Applications for Assistantships, requests for copies of the pub- 
lications and inquiries for information on engineering and in- 
dustrial research problems should be addressed to the Director, 
Engineering Experiment Station, University of Washington, 
Seattle 5. 
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IN RETROSPECT: HOW THE FIRST ATOMIC BOMB AFFECTED 
ONE OF ITS CREATORS FIFTEEN YEARS AGO* 


PRESENT OVERALL STATUS 


As the result of the labors of the Manhattan District or- 
ganization in Washington and in Tennessee, of the scien- 
tific groups at Berkeley, Chicago, Columbia, Los Alamos, 
and elsewhere, of the industrial groups at Clinton, Han- 
ford, and many other places, the end of June 1945 finds us 
expecting from day to day to hear of the explosion of the 
first atomic bomb devised by man. All the problems are 
believed to have been solved at least well enough to make 
a bomb practicable. A sustained neutron chain reaction 
resulting from nuclear fission has been demonstrated; the 
conditions necessary to cause such a reaction to occur ex- 
plosively have been established and can be achieved; pro- 
duction plants of several different types are in operation, 
building up a stock pile of the explosive material. Al- 
though we do not know when the first explosion will occur 
nor how effective it will be, announcement of its occurrence 
will precede the publication of this report. Even if the first 
attempt is relatively ineffective, there is little doubt that 
later efforts will be highly effective; the devastation from 
a single bomb is expected to be comparable to that of a 
major air raid by usual methods. 

A weapon has been developed that is potentially de- 
structive beyond the wildest nightmares of the imagination; 
a weapon so ideally suited to sudden unannounced attack 
that a country’s major cities might be destroyed overnight 
by an ostensibly friendly power. This weapon has been 
created not by the devilish inspiration of some warped 
genius but by the arduous labor of thousands of normal 
men and women working for the safety of their country. 
Many of the principles that have been used were well 
known to the international scientific world in 1940. To 
develop the necessary industrial processes from these prin- 
ciples has been costly in time, effort, and money, but the 
processes which we selected for serious effort have worked 
and several that we have not chosen could probably be 
made to work. We have an initial advantage in time be- 
cause, so far as we know, other countries have not been 
able to carry out parallel developments during the war 
period. We also have a general advantage in scientific and 
particularly in industrial strength, but such an advantage 
can easily be thrown away. 


PROGNOSTICATION 


As to the future, one may guess that technical develop- 
ments will take place along two lines. From the military 
point of view it is reasonably certain that there will be 
improvements both in the processes of producing fission- 
able material and in its use. It is conceivable that totally 
different methods may be discovered for converting matter 
into energy since it is to be remembered that the energy 
released in uranium fission corresponds to the utilization of 
only about one-tenth of one per cent of its mass. Should 
a scheme be devised for converting to energy even as much 
as a few per cent of the matter of some common material, 
civilization would have the means to commit suicide at will. 

The possible uses of nuclear energy are not all destruc- 
tive, and the second direction in which technical develop- 


* The content of this page is from the conclusion of the first public 
report (printed in 1945) on the development of the atomic bomb. Ex- 
cept for short omissions, this excerpt includes all the “General Sum- 
mary” of H. D. Smyth's report, A General Account of the Develop- 
ment of Methods of Using Atomic Energy for Military Purposes under 
the Auspices United States .S. Government 
Printing Office, 1 


ment can be expected is along the paths of peace. In the 
fall of 1944 General Groves appointed a committee to look 
into these possibilities as well as those of military signifi- 
cance. This [Tolman} committee . . . received a multitude 
of suggestions .. . While there was general agreement 
that a great industry might eventually arise, comparable, 
perhaps, with the electronics industry, there was disagree- 
ment as to how rapidly such an industry would grow; the 
consensus was that the growth would be slow over a period 
of many years. At least there is no immediate prospect of 
running cars with nuclear power or lighting houses with 
radioactive lamps although there is a good probability that 
nuclear power for special purposes could be developed 
within ten years and that plentiful supplies of radioactive 
materials can have a profound effect on scienitfic research 
and perhaps on the treatment of certain diseases in a similar 
period. 
PLANNING FOR THE FUTURE 


During the war the effort has been to achieve the maxi- 
mum military results. It has been apparent for some time 
that some sort of government control and support in the 
field of nuclear energy must continue after the war. Many 
of the men associated with the project have recognized this 


fact and have come forward with various proposals, some . 


of which were considered by the Tolman Committee, al- 
though it was only a temporary advisory committee report- 
ing to General Groves. An interim committee at a high 
level is now engaged in formulating plans for a continuing 
organization. This committee is also discussing matters of 
general policy about which many of the more thoughtful 
men on the project have been deeply concerned since the 
work was begun and especially since success became more 
and more probable. 


THE QUESTIONS BEFORE THE PEOPLE 


We find ourselves with an explosive which is far from 
completely perfected. Yet the future possibilities of such 
explosives are appalling, and their effects on future wars 
and international affairs are of fundamental importance. 
Here is a new tool for mankind, a tool of unimaginable 
destructive power. Its development raises many questions 
that must be answered in the near future. 

Because of the restrictions of military security there has 
been no chance for the Congress or the people to debate 
such questions. They have been seriously considered by all 
concerned and vigorously debated among the scientists, and 
the conclusions reached have been passed along to the high- 
est authorities. These questions are not technical questions; 
they are political and social questions, and the answers 
given to them may affect all mankind for generations. In 
thinking about them the men on the project have been 
thinking as citizens of the United States vitally interested 
in the welfare of the human race. It has been their duty 
and that of the responsible high government officials who 
were informed to look beyond the limits of the present war 
and its weapons to the ultimate implications of these dis- 
coveries. This was a heavy responsibility. In a free coun- 
try like ours, such questions should be debated by the 
people and decisions must be made by the people through 
their representatives. This is one reason for the release of 
this report. It is a semitechnical report which it is hoped 
men of science in this country can use to help their fellow 
citizens in reaching wise decisions. The people of the 
country must be informed if they are to discharge their 
responsibilities wisely. 


At the control and instrument console for the new UW Training Reactor 


HE ENGINEERING EXPERIMENT STATION 
announces the availability of research assistant- 
ships for 1961-62. Appointees must qualify for full 
graduate standing in the University. A graduate 
assistant receives $2,700.00 over four quarters at the 
rate of $225.00 per month; he pays tuition, but receives 
a waiver for the difference between resident and non- 
resident fees. These research assistantships require 
that half time be devoted to research leading to a thesis. 

Research opportunities are offered in all branches 
of engineering, including nuclear engineering. (See 
the University of Washington catalog.) Typical 
examples of current research will be found in this and 
other issues of The Trend in Engineering. Assignments 
are open in connection with government, industry, 
and Engineering Experiment Station research projects. 
Teaching assistantships are also open in some depart- 
ments. 

Additional information and application forms may 
be obtained from the Director of the Engineering 
Experiment Station, University of Washington, Seattle 
5, Washington. Applications may be submitted at any 
time, although most appointments become effective 
September 16, 1961. 


Research Assistantships for 1961-1962 


Aeronautical 
Chemical 
Civil 
Electrical 
Mechanical 
Mineral 


Ceramic 


Metallurgical 


Nuclear 
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Plane Elasticity Problems and the Direct Stiffness Method 


HAROLD C. MARTIN 
Professor of Aeronautical Engineering 


Plane elasticity problems 
involve the determination of 
stresses and deformations in 
a two-dimensional continu- 
um. Loads and deformations 
are restricted to those lying 
within the plane of the con- 
tinuum. 

Solutions usually 
achieved by solving the ap- 
propriate differential equa- 


ary conditions of the prob- 
lem. In certain simple cases, 
exact solutions may be found. More often, approxi- 
mate solutions are the only ones that can be de- 
termined. Examples may be found in any text on 
elasticity. 

The purpose of this article is to call attention to a 
different concept for finding approximate solutions 
to such problems. This alternate procedure replaces 
the continuum with an assemblage of finite elements 
and then employs the concepts of the direct stiffness 
method in obtaining a solution. It thereby becomes 
feasible to analyze very complex problems, and 
therein lies the real advantage of this method. Digi- 
tal computing equipment must be used in carrying 
out the numerical calculations. 

The ideas underlying the method were developed 
as part of a general procedure for analyzing complex 
engineering structures. (See Reference 1, in which 
the theory used in this present article was first 
given.) Details and points of view have been modi- 
fied to some extent in the present writing. 


Direct Stiffness Method 

The differential equations of the continuum are 
useful in reaching a solution only if the problem is 
quite simple. Engineering structures built up of bars, 
beams, sheet elements, and so on are generally much 
too complex to be analyzed in this manner. It then 
becomes necessary to apply the fundamental condi- 
tions of equilibrium and continuity in such a manner 
that the mathematical formulation is given in terms 
of algebraic equations. The direct stiffness method 


H. C. Martin 


offers a convenient scheme for obtaining these equa- 


tions. 
In applying the stiffness method, we assume that 
the built-up structure can be analyzed as an assem- 


blage of structural units connected at certain selected . 


points called nodes. The simplest example is the pin- 
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tions, subject to the bound- 


jointed truss. More general structures will include 
beams, plate elements, and so forth as parts of the 
assemblage. However, even in these cases, the vari- 
ous subelements will be assumed to be connected only 
at the node points. Later the two-dimensional con- 
tinuum will also be treated in this same manner. 

Before the stiffness method can be applied to the 
overall structure, the stiffness of the structural sub- 
units must be known. This can be illustrated by 
reference to the typical truss member. If in Fig. l,a 
L = length of member, A = cross-sectional area, E 
= modulus of elasticity, and X and u are force and 
displacement as shown, elementary theory then 
proves that 


In Eq. (1) equilibrium is satisfied. 

For a general orientation of the member as shown 
in Fig. 1,b, appropriate coordinate transformation 
applied to Eq. (1) leads to 


xX, Symmetric | 
Yi AE | ™ 
@) 
X. 2 Nu U2 


where A = cosine of angle between member 1-2 and 
OX, » = cosine of angle between member 1-2 and OY, 
u, = displacement of node 1 parallel to OX, v, = dis- 
placement of node 2 parallel to OY, and so on. In 
compact matrix form Eq. (2) becomes 


{X} =[K]{u}, (3) 


in which [K] is the stiffness matrix for the member. 

If the single bar of Fig. 1,b is considered as a 
structural entity, three displacements in Eq. (2) will 
have to be specified as zero in order to prevent rigid 
body motion. The corresponding forces then become 
reactions. Only one active force component can be 
applied to this member. However, the reactions and 
the active displacement component can be readily 
found from Eq. (2). 

Now consider the truss consisting of m pin-con- 
nected members. Let the individual stiffnesses of 
these be K™, K®,...K™, Then, for the entire 
truss, 


2 X23 


Fic. 1. AXIAL-FORCE MEMBER 


0O 
0 
{X} = {u} (4) 


K™ 

This last equation can be consolidated. The process 
is not difficult, but need not be explained here. The 
salient point is that the stiffness of the truss, 
|K |ervss, follows directly from the square matrix 
of Eq. (4). Consequently, the stiffness matrix for 
the overall structure can be readily obtained provided 
the stiffnesses of the individual members are known. 
Furthermore, a series of matrix calculations then 
leads to the desired displacements and reactions. 

There remains the problem of determining the 
forces of the internal members. These do, of course, 
become known when the nodal displacements have 
been found. In fact, Eq. (2) will enable, for example, 
X, and Y, to be calculated for member 1-2 in terms 
of u,, %, M., and v,. Then the member force S is 
simply given by 


S = XA+U (5) 


All unknowns thereby become known. Note also that 
it makes no difference what the order of redundancy 
may be, nor do we have to perform a stress analysis 
for unit-applied forces or redundants. 


Two-Dimensional Continuum 

The rectangle of Fig. 2,a is a simple example. 
Forces and support conditions have not been in- 
cluded. 

Figure 2,b shows the same rectangle replaced by a 
total of 4 triangular subelements joined together at 
nodes 1,2, .. . 6. No connection between adjacent 
triangles is assumed except at the nodal points ; and 
the selection of triangles is in no sense unique. 
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Fic. 2. Two-DIMENSION CONTINUUM 


When this idealization is kept in mind, the con- 
tinuous rectangle becomes quite similar to the previ- 
ous truss. If we can determine the stiffness matrix 
for an arbitrary triangular unit, we can then form 
[K] for the idealized surface of Fig. 2,b. A complete 
solution is then possible. This will be an approxima- 
tion to the true solution for the actual problem. As 
more triangles are used to represent the given sur- 
face, the accuracy of the results will improve. Ex- 
perience indicates that satisfactory data can be ob- 
tained with a reasonable number of subunits. 

Many additional statements can, and probably 
should, be made. For example, others besides tri- 
angular shapes can be used for the subunits, although 
triangles are the simplest to handle and give accept- 
able results. The density of nodes can also be varied 
so as to be a maximum in regions of maximum com- 
plexity or interest to the analyst. Great flexibility is 
therefore available when replacing the continuous 
surface with a finite assemblage of subunits. 


Stiffness Matrix for the Triangle 

The problem is to find the stiffness matrix [K], 
which relates the six force components Y,, Y,,... to 
the six corresponding displacements, “,, v,, .. 
This is an interesting problem, for which the solution 
first appeared in Reference 1. The arbitrary triangu- 
lar element is illustrated in Fig. 3. 

It is evident that the concentrated node forces will 
cause stress singularities in the triangle. These will, 
of course, be at variance with the usual state of stress 
existing in the triangular subelement replacing the 
continuum. Therefore the nodal forces must be 
viewed as a useful but fictitious representation for a 
reasonable state of stress in the element. 

To accémplish this we begin by assuming constant 
strain in the triangle, or 


and 
Yan? =. (6) 
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Fic. 3. ARBITRARY TRIANGULAR ELEMENT 


Integration gives, 
u=ax+(c—d)y+f 

and (7) 
v=dx+by+e. 


Constants d, e, and f are due to integration. The six 
constants agree with the total number of nodal dis- 
placements. This is necessary to the analysis that 
follows. 


We note that when e,=€,=72y=0 
u=—dytf, 
v= dx-+e, 
and 


ao ou 
Ox dy 


ll 


2a = 


a=b=c=0 


Hence d represents rigid body rotation, while e and f 
represent rigid body translation in the absence of 
rotation. 


It is now possible to write u, v at the nodes. This 
gives six equations that can be solved for the six 
constants in terms of nodal displacements and co- 
ordinates locating the vertices of the triangle. Fur- 
thermore, the strains can now also be expressed in 
terms of nodal displacements. Let this equation be 
given by 


{e} =[@] {x}. (8) 


Hooke’s law relates stresses and strains. The equa- 
tions will be expressed as 


{o} (9) 
In the above equations only constant elements 
appear in [@] and [¢]. 


The strain energy for the triangle can now be 
written as 


where the integration is over the volume of the ele- 
ment. Substituting into this last equation from Eq. 
(8) and Eq. (9) gives 


Um [EV [@] {u}d(vol). 


The integrand is a constant ; hence U is known. 


(10) 


Applying Costigliano’s Theorem permits nodal 
forces per unit displacement to be calculated. This 
determines the required stiffness matrix as follows: 


| Lol let’ [0] x Vol triangle. 
{u} = {1} 


As given in Reference 1, the result of the calcula- 
tion just completed is the stiffness matrix for the 
arbitrary triangle of Fig. 3. As calculated from Eq. 
11, this matrix is as follows: 


(11) 


uy V1 Ve Ug V3 
Et 

[K]= 7? X13" +413 — 
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i 
| 


1 
where 
A =area of triangle 
t =thickness of triangle 
1 


Xig =X;—x,, and so on. 


Since |K] is symmetric, the elements represented by 
can be written at once by inspection. 

The matrices required in stress calculations as 
given by Eq. 12 are 


v 0 
E 1 0 
1-? 
0 0 
= 
and, 
0 0 Vie 0 
1 
0 Xx x13 O 


where, as before, 
and 


Once again we need to keep in mind that the nodal 
forces associated with the arbitrary triangle are in a 
sense fictitious. This is true to the extent that, acting 
through the corresponding nodal displacements, they 
produce a state of constant strain (and stress) in the 
triangle. 

An important point, which was not mentioned in 
Reference 1, is that, since « and v are linear in x and 
y, Eq. (7), common edges of adjacent triangles will 
remain common after deformation has occurred. 
Hence continuity is preserved between nodes. How- 
ever, strains are restricted as given by Eq. (6). 

Finally, having found [K] for the assemblage of 
triangles to determine node displacements, the prob- 
lem remains to find the stresses. This is readily 
accomplished by substituting Eq. (8) into Eq. (9), 
to get 


{o} =[Ello] {x}, (12) 
where [€] and [@] are known. 


The stresses thereby determined for any triangular 


Db 3 
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Fic. 4. TETRAHEDRON 


subelement are uniform over the triangle. However, 
stresses will vary between triangles. Various schemes 
will occur to the analyst for using these results for 
plotting or otherwise estimating the continuous stress 
variation in the actual member. 


Three-Dimensional Continuum 


Before proceeding to examples illustrating the use 
of the stiffness matrix for the triangular element, 
extension of the previous theory to the solid will be 
briefly discussed. Although these extensions have 
not been worked out in detail as yet, they are logical 
developments beyond the planar case. 

The basic element is taken as the tetrahedron (Fig. 
4). A total of 12 components of nodal displacement 
can now exist. If we again choose constant strains, 
we have 


°° ™ Ox 
ov Ow , Ov 
and 
ow Ow , Ou 


Integrating to obtain displacements leads to 
v=(d—g)x+by+(e—h)z+m , (14) 

and 
w=kx+hy+cz+n. 

The above displacement equations contain 12 


arbitrary constants. Consequently the chosen strain 
expressions will enable the stiffness matrix to be 


(Continued on Page 19) 
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A Rational Analytic Theory of Fatigue 


PAUuL C. Paris 
Assistant Professor of Civil Engineering 


Mario P. GOMEz* and WILLIAM E. ANDERSON 
Research Engineers, Boeing Airplane Company 


P. C. Paris 


A great deal of effort has recently centered around 
examination of the factors influencing the growth of 
fatigue cracks. Fatigue has been considered a multi- 
phase problem: e.g., initiation of a crack and its 
growth are often considered as separate phenomena. 
In contrast, the objective of this work is to show that 
the growth of an initial “crack-like” imperfection to 
a critical size, which causes static failure of a struc- 
ture, may be described by a single rational theory. 

Two loading parameters, the nature of the stress 
field near the tip of a crack and the variation of this 
field, are taken to control the rate of crack extension 
in a given material. This hypothesis is proven by 
using it to correlate data from three independent in- 
vestigators. Since it shows a positive correlation of 
all available data for crack-extension rates from 10-7 
to 10°? in. per cycle, the hypothesis may be used to 
formulate a theory of fatigue that permits computing 
the structural lives of complicated geometries from 
simple laboratory tests of material properties. 


The Stress Distribution Near the Tip of a Crack 

The form of the stress distribution in the vicinity 
of a crack root was given by Sneddon’ in 1946 and 
has recently been expanded by Irwin®* and Wil- 
liams.* The unique character of this form, as Irwin 
showed,’ is a controlling factor in attempts to analyze 
crack extension under static loads. We will show that 
this same character becomes fundamental in crack 
extension under cyclic loading upon the addition of 
new concepts to describe the cyclic nature of the 
loading. 

* Mr. Gomez received his M. S. degree in Metallurgical 
Engineering in 1958 at the University, after which he worked 


for Boeing. He is now Senior Scientist at the Missile Sys- 
tems Division of Lockheed Aircraft Corporation. 
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M. P. Gomez 


W. E. Anderson 


Restricting this discussion to cracked bodies in 
which the geometry and loading of the body are sym- 
metric with respect to the plane of the crack results 
in very little loss of generality. The nature of cracks 


y 


Fic. 1. COORDINATES USED To DESCRIBE STRESSES 
NEAR A CRACK TIP (9; ¢,; 7,; v) 


is to form most often on such planes, i.e., planes per- 
pendicular to maximum-principle tension stresses. 
Williams* and Irwin® have given the required forms 
of stresses for other cases, but these will not be dis- 
cussed further in this work. 

The coordinates of points in a cracked body with 


9 
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SYMBOLS 

a =the half crack length 

b = plate dimension parallel to a crack 

F =a loading force on a body 

f() or FQ) =a function of 

K =the stress singularity-intensity factor 
at the tip of a crack 

| = =the critical value of K for a material 
associated with a crack extension 
under static load 

L = the length of a plate 

N =the number of load cycles since 
initial loading 

P = the loading of a body 

r,0 = polar coordinates from the crack tip 

x,y = rectangular coordinates centered 
with respect to a crack 

x =the correction factor for K in plates 
of finite width 

8 = the ratio of maximum to minimum 
load on a body during a load cycle 

= Poisson's ratio 

Sy = gross area stress or nominal stress level 

= components of stress near a crack tip 

a(x) =the normal stress present at a crack 


location before the crack appeared 
Aa/AN or da/dN =the rate of crack extension 


crack-plane symmetric geometry may be described as 
in Fig. 1. If terms of higher order in r are ignored, 
the elastic solution for stresses in the vicinity of the 
crack tip for all such problems is 


= 1 sin sin ] 


(1) 


ta 0, 


and 


ll 


o, = 0(plane stress) 
or 


o. (plane strain) 
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Fic. 2. A TYPICAL FATIGUE-PANEL CONFIGURATION 
FOR WHICH TEST DATA WERE AVAILABLE 


This result implies that the distribution of elastic 
stress always has the same functional form near the 
singularity caused by the sharp crack root and that 
it differs only by a stress singularity-intensity factor, 
K. This factor is linearly dependent upon the loads 
on the body and also must contain a geometric factor 
related to the crack length and other geometric prop- 
erties of the body. An example is K = oV\/a for the 
problem of a crack of length 2a in an infinite sheet 
under a uniform tensile-stress field, o, perpendicular 
to the crack. 


Now the Griffith-Irwin theory of static strength of 
bodies containing cracks may be resolved from this 
discussion in the following fashion: Identical inten- 
sity factors of elastic stress will result in identical 
yield zones. near the tips of cracks in the same ma- 
terial if the yield zones are small compared to the 
region of applicability of the stresses given by Eq. 
(1). As discussed in previous works,® the size of the 
yield zone may be shown to be small if the nominal 
stresses in the body are well below the yield point. 
Therefore, regardless of the appearance of a small 
yield zone, there will be some critical value of the 
stress singularity intensity, K,,, near the tip of a 
crack that will cause static crack extension in a given 
material. The above hypothesis is the equivalent of 
the Griffith-Irwin theory, which was originally based 
on energy considerations. 
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The Significance of Stress-Intensity Factors 
in Fatigue 

The stress-intensity factor may be considered to 
be a measure of the effect of the loading and the 
geometry of a body on the stress intensity near the 
root of a crack. Therefore, as the loads on a body 
vary and as the geometry changes by crack exten- 
sion, the instantaneous values of K reflect the effects 
of these changes at the crack root. 

Let B be the ratio of maximum to minimum load 
on a cracked body during a cycle of loading. Then, 
since K is directly proportional to the magnitude of 
the load, B is also the ratio of Kmax to Kmin, regard- 
less of the geometry of the body; that is, 


Therefore the stresses near the root of a crack are 
completely described by Kmax and B in a given ma- 
terial, since these two parameters give both the in- 
tensity and variation of the effects of loading and 
geometry. 

A theory of fatigue crack extension may now be 
hypothesized as follows: Since, as has been shown, 
during a cycle of loading the stresses and strains near 
the tip of a crack are completely specified by Kmax 
and B, we can reasonably assume that any phenom- 
ena occurring in this region are controlled by these 
parameters. The amount of crack extension per cycle 
of loading is just such a phenomenon, or, in func- 
tional form, 
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Experimental Evidence 


It is pertinent to examine the above hypothesis in 
the light of available experimental data. Many in- 
vestigators’***!° have measured crack-extension rates 
due to cyclic loading in aluminum alloys by using the 
configuration shown in Fig. 2. The majority of these 
tests have been performed with minimum loads near 
zero, or B = «; therefore from examination of Eq. 
(3), these results for a given material should form 
a single curve on a plot of Kmax vs Aa/AN. For such 
a configuration Kyax may be computed from* 


K = 


_ V4+2(a/b) 
where (a/b)? — (a/b»* (4) 
if a<b. 


The results of this attempt at correlation of crack 
extension-rate data are shown on Figs. 3 and 4. It is 


worthy of special note that the data on these curves 
are from three independent investigators, using many 
specimen sizes, i.e., widths from 1.8 to 12 in., thick- 
nesses from 0.032 to 0.102 in., and lengths from 5 to 
35 in. The testing frequencies varied from 50 to 2000 
cpm, and the maximum stresses on the gross area 
varied from 6 to 30 ksi. On each graph, the materials 
are both clad metals and bare metals. Therefore the 
correlation shown is surely more than coincidental. 

On the presumption that such curves may be ob- 
tained for various values of B for a given material 
from laboratory tests, this discussion proceeds to 
formulate the necessary elements of an analytic the- 
ory of fatigue, these results being applied in the fol- 
lowing sections. 


An Analytic Theory of Fatigue 


Knowledge of the material-property curves in the 
form of Figs. 3 and 4, with the addition of curves for 
other B values, implies the functional form of Eq. (3) 
as given. Further, given the loading and geometry 
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of a structure, B is known from the load ratio during 
any cycle, N, or 


6 = 6 (5) 


Moreover, K may be computed for any crack length, 
a, and the maximum load as given during the Nth 
cycle, or 


= K max (N,a). (6) 


Therefore, for a problem with the above specifica- 
tion, Eq. (5) and Eq. (6) may be substituted into 
Eq. (3) to give 


da/dN = F(N,a), (7) 


where the functional form of F is known, point by 
point, from the data given, and effects of loading 
history are neglected. 


The solution to Eq. (7) may be found provided an 
initial crack size may be specified, or some equivalent 
condition may be stated. In practice, maximum im- 
perfection sizes may be stipulated on the basis of 
material quality, production methods, and inspection 
technique, for in practice we know that fatigue cracks 
grow from just such imperfections. Then Eq. (7) 
may be integrated, at least by numerical procedures, 
to generate a complete crack history for the structure. 
The most difficult phase of this analysis is the com- 
putation of K for a given load and crack length, as 
will be commented upon later. 


Relationships to Classical Fatigue Theory 


Suppose the material-property curves of the form 
of Figs. 3 and 4 are known for a specimen that has 
been subjected to an ordinary fatigue test and has 
developed an easily measurable crack in a given num- 
ber of cycles. By using these curves and the measured 
crack length, the crack-extension rates may be inte- 
grated backward to determine an effective initial im- 
perfection size. 

Using this computed initial imperfection size makes 
it possible to calculate the number of load cycles re- 
quired for failure of the specimen at any stress level. 
Therefore the complete S-N curve for a material may 
be computed by this process; moreover, the whole 
S-N curve can be obtained from a single specimen, 
since curves of the material properties for several B 
values may be generated during the same test. 

An accumulative damage theory is automatically 
present in the preceding analysis, which replaces 
Miner’s empirical hypothesis of damage. The Soder- 
berg or Goodman diagrams being likewise empirical, 
this paper presents a rational analytic theory. As is 
evident, the form of the analysis lends itself well to 
statistical analyses. 
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Computation of K, the Stress-Intensity Factor 


By use of data on materials, in the form of Fig. 3, 
the inherently nonlinear problem of structural life 
has been resolved to require only the computation of 
instantaneous values of K as the crack propagates 
through the structure. Since K is the elastic-stress 
intensity factor for stresses near a crack tip, the prob- 
lem has been essentially linearized from the point of 
view of mechanics, a method that permits the use of 
well-known techniques in attacking the problem. 

Irwin*® has given the solution to the problem shown 
in Fig. 5, which may be applied by superposition tech- 
nique (as a Green’s Function) to solve all problems 
of sheet-skins containing internal cracks. This meth- 
od, which is described in detail in Reference 6, will 
be outlined here. 

Consider the general and arbitrary x-axis symmet- 
ric plane-stress problem containing a crack, illus- 
trated by Fig. 6(a). This may be considered to be 
a free body removed from a gross structure, and may 
be chosen large enough in size that the crack itself 
has very little influence on the magnitude and distri- 
bution of the boundary forces and stresses shown. 
Thus we may proceed to solve the plane-stress prob- 
lem for this sheet, but we find considerable difficulty 
in problems of doubly connected regions. As an al- 
ternative, this problem may be considered to be the 
sum of two problems, i.e., those in Fig. 6 (b) and (c). 
The first, (b), is the solution to the same problem 
with no crack present; the second, (c), is the solu- 
tion with the same geometry as the original problem, 
(a), but with loads on the crack surfaces only, of 
equal and opposite intensity to the stresses that occur 
at the crack location in (b). This approach ensures 
that the sum of the boundary forces in (b) and (c) 
are identical to those in (a), and thus their sum is the 
solution to (a). 

Now we desire only to determine the stress inten- 
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Fic. 6. THE REDUCTION OF A PROBLEM, (a), INTO Two 
SIMPLER PROBLEMS, (b) and (c), FOR COMPUTATIONS OF 
STRESS SINGULARITY - INTENSITY FACTORS 


sity factor, K, for the problem, which is the sum of 
the stress intensity factors, K’ and K”, for problems 
(b) and (c). Since K is the intensity of the singu- 
larity of stresses at the crack tip—note 1/\/r in 
Eq. (1)—and this singularity is not present, in (b), 
then 


kK’ =0. (8) 


The task remaining is to determine K” for problem 
(c), since, observing Eq. (2), 


K = K". (9) 


Thus the solution to all problems of the type of (a) 
is reduced to finding the stress, 0, (x), on the crack 
axis with no crack present, as in (b), and using this 
stress as the loading in (c). 

To find the value of the stress singularity, K”, for 
(c), Irwin’s solution* for the stress singularity in 
Fig. 5 gives 
P(a+b)! 


Kp = 


14 


Taking (10) as a Green’s Function, the stress sin- 
gularity, K”, is given by 


which may be further simplified by y-axis symmetry, 


Hence all the sai problems of the form of Fig. 
6(a) may be attacked by solving the problem for the 
stresses, o, (+), along the crack axis with the crack 
absent, and integrating the result by Eq. (11) or 
Eq. 12. 

Equally powerful techniques may be presented for 
other general classes of problems. Therefore we con- 
clude that, in general, K may be computed easily in 
most of the problems of interest for engineering 


purposes. 


Conclusion 


On the basis of the experimental data given, it is 
evident that rates of crack growth—for example, 
those in 2024-T3 and 7075-T6 skins of aircraft struc- 
ture—may be computed by the theory presented over 
a wide range of nominal stress levels and crack sizes. 
The ramifications of such broad correlation imply an 
analytic theory of fatigue based on a concept of 
growth from initial imperfections through which 
structural life may be predicted. 
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High-Temperature Metallography 


B. KASEN 
Research Assistant, Metallurgical Engineering 


Metallography the 
branch of metallurgy that 
relates the microstructure 
of metals and alloys to their 
physical and mechanical 
properties. Conventional 
metallographic techniques 
have provided the metallur- 
gist with much valuable in- 
formation concerning solid- 
state reactions that accom- 
pany physical and mechani- 
cal property changes in ma- 
terials subjected to unusual 
environments. However, with conventional tech- 
niques it has been necessary to base predictions of 
high-temperature behavior on metallographic obser- 
vations made at room temperature. The development 
of practical methods for direct examination of the 
microstructure of metals at elevated temperatures 
has provided a new approach to the study of the 
high-temperature constitution of alloys and of micro- 
structural processes. 

The structure observed under the metallurgical mi- 
croscope must be truly representative of the parent 
material if the conclusions drawn from the study are 
to be valid. Traditionally, the procedure involves 
preparation of an optically smooth surface followed 
by selective etching to reveal structural details such 
as crystal boundaries and phase distributions. Since 
the properties of many alloys are dependent on struc- 
tural changes that occur at elevated temperatures, 
studying the structures that exist in those tempera- 
ture ranges would be desirable. In such studies, com- 
mon practice is to heat a specimen to the temperature 
range desired, hold it until equilibrium is essentially 
achieved, and then rapidly quench the specimen to 
room temperature. The assumption is made that the 
rapid rate of cooling causes the high-temperature 
structure to be retained at room temperature, thus 
allowing the examination to be accomplished. 

The shortcomings of such a procedure are evident. 
The microstructure observed at room temperature 
after quenching may not represent the high-tempera- 
ture structure if a phase transformation occurs dur- 
ing cooling. Martensitic transformations or rapid 
precipitation reactions may completely mask the ele- 
vated temperature constitution of an alloy. In addi- 
tion, the preparation of a metallographically smooth 
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surface usually involves 
mechanical polishing, which 
can alter the structure of 
some metastably retained 
phases by deforming the 
surface layers. 

Owing to the problems 
outlined above, it is under- 
standable that metallogra- 
phers have long wished for 
a facility that would permit 
direct observation of the 
microstructure at elevated 
temperatures. In addition 
to eliminating the sources of error just mentioned, 
such a facility would permit continuous observation 
of the same field of view over a wide range of tem- 
peratures. Thus the mechanisms underlying struc- 
tural changes could be studied by direct observation 
of the microstructural changes accompanying phase 
transformations. Such a facility, called a “hot-stage” 
metallograph, is now in operation in the Division of 
Metallurgical Engineering of the School of Mineral 
Engineering. The equipment is currently being used 
in a study of phase transformations in the copper- 
silicon alloy system. This paper will describe the 
facility and illustrate the use presently being made 
of it. 


D. H. Polonis 


Principles of High-Temperature Metallography 


Before proceeding with a description of the high- 
temperature metallograph, we should review some 
basic metallographic principles. It will be subse- 
quently shown that these considerations are influen- 
tial in the design of high-temperature metallographic 
equipment. 

When ordinary “bright field” illumination is used, 
the light beam strikes the specimen essentially nor- 
mal to the surface, and observation is made by means 
of light reflected into the objective. When the sur- 
face has elevations and depressions resulting from 
selective attack by etchants, these structural details 
can be resolved as a result of differences of intensity 
of the light reflected from the surface irregularities. 
For example, the presence of a surface elevation may 
indicate the presence of a phase not greatly attacked 
by etchants, whereas a valley may be due to a phase 
with a higher rate of solution or to a crystal bound- 
ary. Thus, for bright-field examination, the surface 
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contour must be altered so that differences in the in- 
tensity of light reflected from the structural details of 
the surface will enable the observer to resolve these 
details and to see the microstructure. 

If the light incident on a metallic surface is plane 
polarized and if the material under study is isotropic 
(for example, a phase of cubic structure), then the 
reflected light will also be plane polarized. By proper 
placement of a polarizing filter in the microscope ob- 
jective, extinction of reflection from such surfaces can 
be obtained. However, if the structure also contains 
an anisotropic phase (for example, a phase of hex- 
agonal structure) then complete extinction will not 
occur because light reflected from such phases is no 
longer plane polarized. A resolution of separate 
phases may thus be achieved. 

The principles by which a structure is revealed 
under bright field or under polarized light are appli- 
cable regardless of the temperature at which the ob- 
servation is made. If bright field illumination is used, 
the surface must have a contour that reflects the un- 
derlying structure or, if polarized light is used, the 
desired distinction must be between isotropic and 
anisotropic phases. 

There are essentially three ways in which heating 
can produce contour on a metallographically smooth 
metal surface: (1) thermal etching, (2) surface tilt- 
ing, and (3) surface upheaval.* 

Thermal etching occurs when a metal or alloy is 
heated to a sufficiently high temperature in a vacuum 
or in a protective atmosphere. It occurs where a 
grain boundary or an interface between phases inter- 
sects the surface of the specimen. It can also occur 
where a plane of low surface energy exists near the 
surface. The mechanism of grain-boundary etching is 
believed to involve the establishment of equilibrium 
surface tension between the intersecting boundaries. 
Thus the angles alpha and beta of Fig. la are estab- 
lished by means of surface diffusion or by evapora- 
tion from the boundary and deposition adjacent to it. 
In a somewhat similar fashion, net evaporation may 
cause low-surface-energy planes to develop by trans- 
port of metal through the vapor phase, provided that 
means exist for removal of the vapor from the sys- 
tem.? This may be accomplished by a dynamic vac- 
uum system or by condensation on cooler parts of a 
furnace. 

Surface tilting and surface upheaval are similar in 
that both reflect a change in volume, concomitant 
with the phase change. The distinction is made on 
the basis of the shape of the surface after the trans- 


formation. Surface tilting, illustrated in Fig. 1b, oc- . 


curs when the transformation is of the “martensitic” 
type. In this transformation, rows of atoms suddenly 
shear across each other, producing a rotation in addi- 
tion to a volume increase. Surface upheaval, illus- 


16 


a. THERMAL ETCHING 


b. SURFACE TILTING 


SURFACE UPHEAVAL 


Fic. 1. THE THREE WAYS IN WHICH HEATING PRODUCES 
CONTOURS ON SMOOTH METAL SURFACES 


trated in Fig. lc, reflects a transformation known as 
a “nucleation and growth” type. The surface contour 
produced by this type of transformation is not well 
defined because of variations in types of transforma- 
tion products, in their shapes and in the locations at 
which nucleation occurs. 

The microstructure may thus be revealed without 
the use of chemical etchants, provided that either 
thermal etching takes place or that a volume change 
occurs. While a specimen cannot be etched under 
vacuum at elevated temperature by using the ordi- 
nary chemical etchants, sometimes a gas at a low 
pressure can be introduced to selectively attack the 
surface of the specimen and reveal its structure. For 
example, a small amount of oxygen may be intro- 
duced to form oxide layers of varying thickness on 
steel alloys, thus permitting phase resolution. 


Hot-Stage Metallographic Equipment 

The type of equipment required for high-tempera- 
ture metallography is shown in Fig. 2. The basic 
components of the facility include a metallograph, a 
hot-stage furnace, a vacuum system, an inert gas 
system, and a control panel.* 

The metallograph, A, is an inverted microscope of 
the type conventionally used in metallographic 
studies. An eyepiece is provided for visual examina- 
tion and provisions are made for photomicrography 
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Fic. 2. Hot-STAGE METALLOGRAPHIC EQUIPMENT 


with 34x 4% in. film. The metallograph is easily 
adapted for 35-mm still photography or for 16-mm 
motion-picture photography. The specimen may be 
illuminated by a tungsten light source attached to the 
metallograph or by a carbon-arc illumination source 
when higher intensity is required. The arc illumina- 
tor is situated on a track, enabling it to be moved out 
of the way when the tungsten source is in use. 

The hot-stage furnace, B, is bolted to the top of 
the metallograph stage. Its design permits the speci- 
men to be observed through a quartz window in the 
bottom of the stage. We shall subsequently return 
to a more detailed consideration of this unit. 

The vacuum unit, C, consists of a conventional 
diffusional pump and a mechanical forepump. The 
pressure inside the stage must be reduced to a level 
where the rate of oxidation of the specimen is negli- 
gible over the period of time required for the study. 
Ultimate pressures between 10-* and 10-° mm of mer- 
cury are adequate for most studies. This system is 
connected to the furnace through a cold trap and a 
flexible metal hose. 

An inert gas such as argon or helium may be used 
to accelerate the rate of cooling of a specimen that has 
been heated under vacuum, or an inert gas may be 
used to protect the specimen during the entire run. 
When inert gas is used, the pressure is maintained 
at approximately one atmosphere by use of a man- 
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Fic. 3. FURNACE FOR Hot-STAGE METALLOGRAPHY 


ometer-bubbler, D. The inert gas is purified by pass- 
ing it over heated calcium chips before it enters the 
furnace. 

The control unit, E, was designed and built in the 
Division of Metallurgical Engineering. It consists of 
three major sections. The top compartment contains 
a strip-chart recorder on which the temperature of 
the surface of the specimen is both indicated and re- 
corded. The center compartment contains the flow- 
meters for controlling the flow of coolant to the fur- 
nace and to the diffusion pump. A thermocouple 
vacuum gauge and a discharge vacuum gauge are also 
located in this compartment. The first of these is 
used for rough vacuum measurement; the second is 
used to measure the lower pressures. The bottom 
compartment contains the power supply, switches, 
and meters. 

A closer view of the hot-stage furnace is provided 
in Fig. 3. Specimen A, a cylinder 10 mm in diameter 
by 25 mm long, is placed on a quartz window with 
the surface to be studied positioned downward. Pre- 
cise temperature measurement is provided by a plati- 
num/platinum-13% rhodium thermocouple, B, rest- 
ing in a slot in the specimen surface. The specimen 
is heated by a tungsten resistance element, C, which 
is lowered over the specimen when the upper and 
lower portions of the furnace are joined. When op- 
erated under vacuum, the element is capable of heat- 
ing the specimen to 1500° C. Thermal conductivity 
somewhat reduces the maximum obtainable tempera- 
ture when an inert gas is used in place of a vacuum. 

When a metal or an alloy is heated to high tem- 
perature in a vacuum, a temperature-pressure rela- 
tionship may be reached wherein volatilization of the 
specimen becomes appreciable. Excessive condensa- 
tion of volatilization products on the quartz window 
will obscure the sample. Consequently, a shutter, D, 
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has been incorporated into the stage. Closing the 
shutter slips an extra piece of quartz between the 
specimen and the microscope objective; opening the 
shutter removes the extra quartz and any condensate 
from the optical path, making it possible to observe 
the specimen. 


Examples of High-Temperature Photomicrography 


An example of photomicrography at elevated tem- 
peratures is shown in Fig. 4. This sequence, photo- 
graphed at 200 diameters, shows the changes that 
occur in the structure of an as-cast copper-5.3% sili- 
con alloy when it is heated to a temperature where 
incipient melting has occurred. An etchant has been 


820 °C 


used to permit observation of the room-temperature 
structure. 

The appearance of the microstructure at room 
temperature is shown in Fig. 4a. Three distinct 
phases may be observed: (1) The matrix, alpha, is a 
solid solution of silicon in copper and has a face- 
centered cubic crystal structure. (2) The dark con- 
stituent has been identified as gamma, which, accord- 
ing to the phase diagram, may have a composition 
ranging from 8.3 to 9% silicon. Previous investiga- 
tors *° have found that the gamma phase has a com- 
plex crystallographic structure of the beta-manganese 
type. The amount of gamma in the microstructure is 
greater than would be expected for a 5.3% silicon 
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Fic. 4. MICROSTRUCTURAL CHANGES DURING HEATING OF A COPPER-5.3% SILICON ALLOY 
Starting Condition As-Cast and Etched (200X) 
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alloy on the basis of the equilibrium diagram. Non- 
equilibrium conditions during freezing have probably 
caused silicon segregation, accounting for the cored 
appearance of Fig. 4a. (3) The pools surrounding 
the gamma consist of kappa phase, which has been 
identified positively by polarized-light examination 
and by hardness measurements. This phase, which 
has a hexagonal close-packed structure, occurs in the 
copper-silicon phase diagram between 5.2 and 6.9% 
silicon, depending on the temperature. The relatively 
large quantity of kappa also reflects segregation in 
the cast alloy. 

At 700° C (Fig. 4b), some fine precipitate particles 
that are not present at room temperature can be seen 
in the alpha solid solution. The gamma and kappa 
phases are still in evidence. The precipitate particles 
exhibited optical anisotropy and are probably kappa 
phase. At 820° C, (Fig. 4c), a distinct change has 
occurred in the appearance of the microstructure, re- 
flecting the start of “incipient fusion,” i.e., local melt- 
ing of regions enriched in silicon content due to the 
coring effects mentioned previously. The grain- 
boundary fusion is clearly revealed in Fig. 4d. 


Conclusions 


Research facilities in the Division of Metallurgical 
Engineering have been substantially increased 
through the acquisition of the equipment described, 
enabling the microstructure of metals and alloys to 
be examined at elevated temperatures. This equip- 
ment is currently being used in a study of elevated 
temperature transformations in the copper-silicon 
system. Subsequent studies will involve other alloy 
systems. The technique will be extremely valuable 
in advancing basic understanding of the mechanisms 
of phase changes in a variety of alloys. 
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DIRECT STIFFNESS METHOD 
(Continued from Page 8) 


calculated along lines similar to these employed for 
the triangle. Moreover, continuity between common 
faces of adjacent tetrahedra will be preserved after 
deformation due to the linear form of Eq. (14). 

Modern computing equipment has already reached 
a state of development that makes it reasonable to 
employ an assemblage of tetrahedra for analyzing 
three-dimensional problems. 


Examples 


To date, applications of the stiffness method have 
been principally directed toward large-scale aircraft 
structural problems. As a result, problems of the 
continuum, as such, have received only minor atten- 
tion. Nevertheless, a useful paper on the subject has 
already been published. (See Reference 2.) 

Figure 5 illustrates three problems investigated in 
Reference 2. Cases a and b are rather obvious; c 
represents the cross section of a dam. For each of 
these problems, deformations as well as stresses were 
determined by using the stiffness method. 


4 HYDROSTATIC PRESSURE 
- + 
j = 
a b c 


Fic. 5. Two-DIMENSIONAL CONTINUA PROBLEMS 
(See Ref. 2) 


Results for problem a are compared against a 
known series solution. Agreement is shown to be 
quite satisfactory, even when relatively few sub- 
elements were used to represent the continuum. Case 
b illustrates a problem difficult to treat by the meth- 
ods of elasticity. Use of the stiffness method would 
make entirely possible the inclusion of additional com- 
plicating features such as taper in width and thick- 
ness. Problem c is analyzed for both dead loads and 
hydrostatic pressure. The results are shown ‘to agree 
satisfactorily with previously calculated and published 
data based on the relaxation procedure. (See Refer- 
ence 2.) 

Additional calculations, presently unpublished, con- 
firm the conclusions reached in Reference 2, namely, 
that the stiffness method provides a useful and versa- 
tile tool for investigating planar continua problems. 
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, “Further Studies on the Local Deformations of 
Honeycomb Face Plate Subjected to Uniaxial Compres- 
sion, Jour. of the Aero/Space Sciences. To be published. 

———., “On the Local Stability of Honeycomb Face Plate 
Subjected to Uniaxial Compression,” Jour. of the Aero/ 
Space Sciences, Vol. 26, No. 10 (Oct., 1 : 

Mills, B. D., “The Fluid Column,” Amer. Jour. of Physics, 
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(with D. A. Bolstad), The Trend in Engineering, Vol. 12, 
No. 3 (July, 1960). 


M.S. in Ceramic Engineering 


Corwin, R. E., Compositional Modifications of Tungsten 
Trioxide and Their Dielectric Properties (1959). 
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of the submicroscopic behavior of ferroelectric crys- 
tals. A potential byproduct of the work is improved 
tuning for the Maser amplifier, an improvement that 
could lead to more effective tuning of radio signals 
from outer space. Graduate students W. R. Elliott 
of Bellevue, Washington, and Attilio Giarola of 
Brazil will work with Dr. Bjorkstam on this project. 


THE TREND IN ENGINEERING 


| 


PUBLICATIONS 

of the 

ENGINEERING EXPERIMENT STATION 
UNIVERSITY OF WASHINGTON 


Requests for copies of publications should be addressed to the Director, Engineering Experiment Station, University of Washington, Seattle 5. 


48. 
54. 
55. 
56. 


58. 


60. 
61. 


67. 
72. 


74. 


76. 
77. 
81. 
85. 
87. 
88. 


91. 


92. 
94. 
98. 
100. 


101. 


102. 
103. 
104. 
106. 
108. 
109. 
110. 
111. 
112. 


BULLETINS 


The Influence of Pipe Thickness on Re-Entrant Intake Losses. 
Charles William Harris. 33 pp. 19 

Elimination of Hydraulic Bd Current Loss at Intake. 
William Harris. 30 pp. 

Vacuum-Tube Control Wind-Tunnel Balances. 
S. Eastman. 20 pp. 

Constant Flow Chaseatestatic of the Plane Orifice in Proximity to 
Side Walls. Charles William Harris. 18 pp. 1931.. 

Preliminary Tests of Thirteen Washington Coals in a Powdered 
Coal Boiler Plant at the University of Washington. Far Sam- 
uel Wilson, Harry Fagan Yancey, and Joseph Daniels. 27 pp. 
1931. 

An Electromagnetic Balance for Force Measurement of Current 
Control. Fred Scoville Eastman. 29 pp. 

Hydraulic Flow Characteristics of a Square-Edged Intake. Charles 
William Harris. 21 pp. 1932. 

The Doherty Washer Cooler, Sybren Ruurd Tymstra. 61 pp. 1932. 
The Derivation of the Mean Temperature Difference Between a 
Saturated Gas and a Cooling Liquid, and also the Derivation of 
= pom Velocity of Such Gas. Sybren Ruurd Tymstra. 19 pp. 
1933. 

The Deformation of Hard Grade Reinforcing Steel During and 
After Construction in a Reinforced Concrete Building. Sergius 
Ivan Sergev. 33 pp. 1933. 

pnd Water Tube Gas Condenser. Sybren Ruurd Tymstra. 24 pp. 
1934. 


Hewitt Wilson. 184 


Charles 
Fred 


Kaolin and China Clay in the Northwest. 
pp. 1934. (Loan only.) 

A Study of Fundamental Relations of the Mathematical Theory 
of Elasticity. Arthur Melvin Winslow. 25 pp. 1934. (Loan only.) 
Influence of Two Secondary Factors in Weir Measurement. 
Charles William Harris. 16 pp. 5 

Discharge of Three Commercial Cippoletti Weirs. Robert Bowman 
Van Horn. 32 pp. 5. 

Open Water Testing of Ship Models. 

20 pp. 1935. 

The Properties of Quartz Sands Washed from Kaolins of the 
Pacific Northwest. Hewitt Wilson and Frank Joseph Zvanut. 42 
pp. 

The Analvsis and Control of Landslides. 
57 pp. 1936. (Loan only.) 

Moment Equations. A Method of Analysis for Continuous Beams 
and Rigid Frames. Alfred Jensen. 88 pp. 

Evaporation and Surface Cooling. George Samuel Wilson and 
William Lyle Dudley. 18 pp. Ss 

Geology and Its Relation to the Occurrence of Oil in Washington. 
Charles Edwin Weaver. 16 pp. 1938. 

Modernization of a Transit System. Factors that epetecion the 
Choice of a Vehicle. George Lisle Hoard. 36 pp. 
Buckling Criteria and Design of Two-Story A slag 
Sergev. 48 pp. 

An Investigation of Prestressed Stays. T. M. Rowlands and F. 
C. Smith. 36 pp. 

Bismuth patos Magnetic Flux Meter. 
48 pp. 1940 

The Telacheonsmster, Telephone Service Measurement in Units 
of Time. Lyall Baker Cochran. 68 pp. 

The Shearing Strength of Cement Mortar. ‘fetes C. Smith 
and Robert Q. Brown. 32 pp. 1941. (Loan only.) 

Silica Sands of Washington. Hewitt _ Kenneth C. Skin- 
ner, and Albert H. Couch. 76 pp. 1942. 

The Number of Teeth in Contact as a Vibration Factor in Invo- 
lute Cut Gears. Sybren Ruurd Tymstra. 16 pp. 

Solutions for Torsional Vibration Frequencies by Reersy Method. 
Arthur Melvin Winslow. 22 pp. 1943. 

Stabilized Earth ~t for Protective Construction. Robert G. 


ennes. 24 pp. 
17 pp. 1944. 


Thomas McKie Rowlands. 


Robert Graham Hennes. 


Sergius I. 


George Sherman Smith. 


Charles W. Harris. 


Hydraulic Models. 


113. 


114. 


116. 


116. 


116. 


116. 


116. 


118. 
119. 


> 


T3. 


T8. 
T9. 


T12. 


The Theoretical Behavior and Design of Initially Curved Struts 
Under an Intermediate Concentric Axial Load. Sergius 
Sergev. 32 pp. 

The Effect of Shearing Forces on the Deflection and Strength of 
Beams. Sergius I. Sergev. 20 pp. 1947. 

Aerodynamic Stability of Suspension Bridges with Special Ref- 
erence to the Tacoma Narrows Bridge. Part I. Investigations 
Prior to October, 1941. F. B. Farquharson. 103 pp. 1949. (Loan 
only.) 

Aerodynamic Stability of Suspension Bridges with Special Refer- 
ence to the Tacoma Narrows Bridge. art II, Mathematical 
Analyses. F. C. Smith and G. S. Vincent. 63 pp. 1950. 
Aerodynamic Stability of my sion Bridges with Special Refer- 
ence to the Tacoma Narrows Bridge. Part III, The Investigation 
of Models of the Original Tacoma Narrows Bridge Under the 
Action of Wind. F. B. Farquharson. 125 pp. 1952. 
Aerodynamic Stability of > ension Bridges with Special Refer- 
ence to the Tacoma Narrows Bridge. Part IV. The Investigation 
of Models of the New Tacoma Narrows Bridge Under the Action 
of Wind. F. B. Farquharson. 104 pp. 1954. 

Aerodynamic Stability of Suspension Bridges with Special Refer- 
ence to the Tacoma Narrows Bridge. Part V, Extended Studies: 
Logarithmic Decrement, Field Damping, Prototype Predictions, 
Four Other Bridges. G. S. Vincent. 104 pp. 1954. 
Investigations in Aeronautics from the Department of Aeronau- 
tical Engineering. Edited by R. M. Rosenberg. 120 pp. 1951. 
The Extremum Principles of the Mathematical Theory of Elas- 
ticity and Their Use in Stress Analysis. (The Walker-Ames lec. 
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$45,900 for basic research on the spinning action of 
electrons. The research is to be a fundamental study 
of the submicroscopic behavior of ferroelectric crys- 
tals. A potential byproduct of the work is improved 
tuning for the Maser amplifier, an improvement that 
could lead to more effective tuning of radio signals 
from outer space. Graduate students W. R. Elliott 
of Bellevue, Washington, and Attilio Giarola of 
Brazil will work with Dr. Bjorkstam on this project. 
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Charles William Harris. 16 pp. 1935. 

Dineres of Three Commercial Cippoletti Weirs. Robert Bowman 
Van Horn. 32 pp. 1935. 

Open Water Testing of Ship Models. Thomas McKie Rowlands. 

20 pp. 1935. 

The Properties of 

pp. 

The por and Control of Landslides. Robert Graham Hennes. 

57 pp. 1936. (Loan only.) 

Moment Equations. A Method of Analysis for Continuous Beams 

and Rigid Frames. Alfred Jensen. 88 pp. 

Evaporation and Surface Cooling. George Samuel Wilson and 

William Lyle Dudley. 18 pp. 1937. 

Geology and Its Relation to the Occurrence of Oil in Washington. 

Charles Edwin Weaver. 16 pp. 1938. 

Modernization of a Transit System. Factors that Peimneine the 

Choice of a Vehicle. George Lisle Hoard. 36 pp. 

Buckling Criteria my Design of Two-Story ll Sergius I. 

Sergev. 48 pp. 1940. 

An Investigation of Prestressed Stays. T. M. Rowlands and F. 

C. Smith. 36 pp. 1941. 

ge a Magnetic Flux Meter. George Sherman Smith. 

48 pp. 1 

The Telechronometer, Telephone Service Measurement in Units 

of Time. Lyall Baker Cochran. 68 pp. 

The Shearing Strength of Cement Mortar. ‘Frederick C. Smith 

and Robert Q. Brown. 32 pp. 1941. (Loan only.) 

Silica Sands of Washington. Hewitt Wilson, Kenneth C. Skin- 

ner, and Albert H. Couch. 76 pp. 1942. 

The Number of Teeth in Contact as a Vibration Factor in Invo- 

lute Cut Gears. Sybren Ruurd Tymstra. 16 pp. 

Solutions for Torsional Vibration Frequencies by Energy Method. 

Arthur Melvin Winslow. 22 pp. 1943. 

Stabilized Earth Blocks for Protective Construction. Robert G. 

Hennes. 24 pp. 1943. 

Hydraulic Models. Charles W. Harris. 17 pp. 1944. 


uartz Sands Washed from Kaolins of the 
ewitt Wilson and Frank Joseph Zvanut. 42 
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116. 


116. 


116. 


116, 


118. 
119. 


to 


T3. 


T8. 
T9. 


T12. 


The Theoretical Behavior and Design of Initially Curved Struts 
Under an Intermediate Concentric Axial d. Sergius I. 
Sergev. 32 pp. 1945. 

The Effect of Shearing Forces on the Deflection and Strength of 
Beams. Sergius I. Sergev. 20 pp. 1947. 

Aerodynamic Stability of Suspension Bridges with Special Ref- 
erence to the Tacoma Narrows Bridge. Part I. Investigations 
to to October, 1941. F. B. Farquharson. 103 pp. 1949. (Loan 
only.) 

Aerodynamic Stability of Suspension Special Refer- 
ence to the Tacoma Narrows Bridge. {athematical 
Analyses. F. C. Smith and G. S. wg 963 4 1950. 
Aerodynamic Stability of Saupe nsion Bridges with Special Refer- 
ence to the Tacoma Narrows Bridge. Part III, The Investigation 
of Models of the Original Tacoma Narrows ’Brid ge Under the 
Action of Wind. F. B. Farquharson. 125 pp. 1952. 


. Aerodynamic Stability of oy nsion Bridges with Special Refer- 


ence to the Tacoma Narrows Bridge. Part IV. The Investigation 
of Models of the New Tacoma Narrows Bridge Under the Action 
of Wind. F. B. Farquharson. 104 pp. 1954. 

Aerodynamic Stability of Suspension Bridges with Special Refer- 
ence to the Tacoma Narrows Bridge. Part V, Extended Studies: 
Logarithmic Decrement, Field Damping, Prototype Predictions, 
Four Other Bridges. G. S. Vincent. 104 pp. 1954. 
Investigations in Aeronautics from the Department of Aeronau- 
tical Engineering. Edited by R. M. Rosenberg. 120 pp. 1951. 
The Extremum Principles of the Mathematical Theory of Elas- 
ticity and Their Use in Stress Analysis. (The Walker-Ames lec. 
tures of the Spring Quarter, 1950.) W. Prager (Brown Univer- 
sity). 30 pp. 1951. 


REPORTS 


. Iron and Steel Manufacture in Washington, Oregon, California, 


and Utah. Joseph Daniels. 69 pp. 1929. (Loan only.) 


. Coal in Washington. Joseph Daniels. 17 pp. 1934. (Loan only.) 
. Water Resources of Washington. Richard Gaines Tyler. 61 pp. 


1938 (Loan only.) 


. A Natural Uranium-Graphite Subcritical Reactor System for 


Education and Research. Part I. Description a and Experimental 
Measurements. A. L. Babb, L. Grimsruc Wilson. Part 
II. IBM 650 Computer Programs for Eitective "Size, Diffusion 
Length, and Material Buckling Calculations. A. L. Babb, L. 
Grimsrud, S. Tashiro. 33 pp. 1959. 


TECHNICAL NOTES 


. The Theory of Fractional Parts Developed for Structural Beam- 


Column Members Preceded by an Elementary Discussion of 
Simple Beam and A ee Design. Sybren Ruurd Tymstra. 32 
pp. 1934. (Loan only.) 


. Gaging Chart for Horizontal Cylindrical Tanks with Flat Heads 


or with Convex Heads. Sybren Ruurd Tymstra. 4 pp. 


. Gaging Chart for Cylindrical Tanks oe Types of 
1 


Heads. Sybren Ruurd Tymstra. 7 pp. 


TREND REPRINTS AVAILABLE 


A New Critical Phase of the Lake Washington Pollution Prob- 
lem. R. O. Sylvester, W. Thomas Edmondson, R. H. Bogan. 
7pp. 1956. 

bd fa Going to Teach Them and How. Editorial Study. 2 pp. 
A Model Study of Some Fish Facilities for the Preposed Mayfield 
Dam. Part I: The Skimmer. E. P. Richey and H. ii. henoweth. 
8 pp. 1957. 

Professional Ethics and Engineering Education. 
Hemenway. 1 p. 1958. 
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